EL 724 020 159 US 



081225.024 5694 



METHOD OF SIMULATION OF PRODUCTION PROCESS OF 
SEMICONDUCTOR DEVICE AND SIMULATOR THEREOF 

CROSS REFERRENCE TO RELATED APPLICATIONS 
5 The subject application is related to subject matter disclosed in the 

Japanese Patent Application 2000-085268 filed in March 24, 2000 in Japan, 
to which the subject application claims priority under the Paris Convention 
and which is incorporation by reference herein. 

1 0 BACKGR OUND O F THE IN VENTIO N 

1. Field of the Invention 

The present invention relates to technologies for the simulation of 
a semiconductor process including oxidation and film formation, and 
particularly, to a method of the simulation of a topography for a 
15 geometrical change of a material region, a simulator for realizing the 
simulation method and a program for accomplishing the simulation 
method. 

2. Description of the Related Art 

20 In a topography simulation for calculating the topography of a 

semiconductor device in an oxidation, deposition and etching step or the 
like, a string model showing the topography of a semiconductor device as 
shown in FIG. 1A is usually used. In the string model, the topography of 
a semiconductor is shown by a string of a line segment Sj connecting a 

25 boundary point Pi and a boundary point P i+ i and each boundary point P is 
moved in conformity with each process model at each small time interval 
to calculate the hourly variation in device topography. 

Also, as shown in FIG. IB, a mesh provided with plurally divided 
triangle elements E; is disposed in the region enclosed by a string to 

30 calculate oxidation and impurity profiles. The boundary point Pj and a 
point Ni+i within the region are called a node N. In the case where the 
mesh exists in the region, all nodes N including not only the boundary 
point but also points within the region are moved at each small time 
interval to calculate the hourly variation in the topography of a 

35 semiconductor device. A topography simulation method using such a 
string model is disclosed in, for example, "VLSI TECHNOLOGY, S.M. 
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Sze, McGraw-Hill, (1988)". However, the prior technologies have the 
problems shown below. 

FIG. 2 shows a flowchart of a conventional method for controlling 
the length of a line segment s. First, in a step S61, the number ni oop of 
5 steps as to what times the calculation of the displacement is looped is set 
to a required number and the present number j of steps is set to zero. 
Next, in a step S62, 1 is added to the present number j of steps. In a step 
S63, the component of all boundary points P as to the displacement is 
calculated, and in a step S64, all boundary points P are moved. In a step 

10 S65, the length of at least one or more line segments s constituting a string 
is regulated. The aforementioned steps S61-S65 are repeated until the 
present number j of steps is equal to ni oop - 

FIG. 3 is a flow chart showing an internal flow of adding boundary 
point P to regulate the length of the line segment s in the step S65 of FIG. 

15 2. First, in a step S71, the length of the line segment s after all boundary 
points P are moved. Next, in a step S72, the identification number i of 
the line segment s is set to 1. In a step S73, whether or not the length ri 
of the line segment Si after the boundary point P is moved is greater than a 
maximum length l max to be specified is judged. If n > l max , a new 

20 boundary point P is added to the line segment si in a step S74. In a step 
S75, 1 is added to the number i. In a step S76, whether or not the number 
i is smaller than the number obtained by adding 1 to a maximum 
identification number (m-1) is judged. If it is judged to be smaller, the 
flow is returned to the step S73 and these steps S73-S75 are repeated until 

25 the identification number i is equal to m. If the identification number i is 
equal to m, the process is finished. Namely, all line segments s were 
judged in order of the number i as to whether or not the length r is greater 
than the maximum length l max . If the length r of the line segment s is 
larger than the maximum length l max , a new boundary point P is added to 

30 the line segment s. 

FIG. 4 is a flow chart showing an internal flow to eliminate a 
boundary point P to thereby regulate the length of the line segment s in the 
step S65 of FIG. 2. First, in a step S81, the length of all line segments s 
after the boundary point P is moved is calculated. Next, in a step S82, 

35 the identification number i of the line segment s is set to 1. In a step S83, 
whether or not the rength n of the line segment si after the boundary point 
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P is moved is smaller than a minimum length l m j n to be specified is judged. 
If r 3 ^ l m jn, in a step S84, one of the boundary points P constituting the 
line segment s is eliminated and a new line segment s'i is produced. In a 
step S85 and a step S86, when the eliminated boundary P is a connecting 
5 point between the line segments si and S2, S2 is designated as s'i. These 
steps are repeated until the identification number i is equal to m-1. 

However, if the length r of the line segment s is set to the 
maximum length l max to be specified or more, a new boundary point is 
always added to the line segment, giving rise to the problem that useless 

10 boundary points are increased. For example, when a groove 59 is formed 
on a silicon substrate 52 as shown in FIG. 5A, a silicon oxide film 61 is 
formed on the surface of the silicon substrate 52 by deposition as shown in 
FIG. 5B. The result of topography simulation as to a region 60 is shown 
FIG. 6. As shown in FIG. 6, the surface string of the substrate 52 is 

15 constituted of line segments arranged at unequal intervals. Each 

topography between boundary points P0 and PI and between boundary 
points P15 and P16 is considered to be able to be estimated from each 
displacement between the boundary points P0 and PI and between the 
boundary points P15 and P16 respectively and no node is arranged 

20 between the boundary points P0 and PI and between the boundary points 
P15 and PI 6, whereby it is expected to shorten calculation time. 
However, because the value of the specified maximum length l ma x is short, 
there is the case where needless boundary points 62 and 64 are prepared 
and therefore the calculation time cannot be shortened as is expected. 

25 On the other hand, it is effective to dispose a new node 63 on 

lengthened line segments sl3 and sl4 having lengths of rl3 and rl4 
respectively to raise the accuracy of topography simulation. However, 
because the specified maximum length l max may be larger, the node 63 is 
formed on the line segment sl4 but no node is produced on the line 

30 segment sl3 and therefore the accuracy of the topography cannot be 
raised. 

In this manner, the value of an absolute maximum length l max 
which shortens time and improves accuracy does not exist, causing such an 
inconvenience that customers must design an optimum value of each of the 
35 maximum length l max and minimum length l min of the line segment s from 
the structure of a semiconductor which must be analyzed by the customers. 
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Specifically, the conventional topography processing methods involve the 
case where exact topography of a device cannot be expressed and the case 
where calculation is made in vain. Therefore, it cannot be said that these 
method are practical to develop semiconductor devices. 

5 

SUMMARY OF THE INVENTION 

The present invention has been conducted in view of the above 
situation and it is an object of the present invention to provide a 
two-dimensional or three-dimensional process simulation method which 
10 can shorten calculation time and improve the accuracy of calculation even 
in a fine region resulting from, for example, the progress of thin 
film-formation and a process simulator which attains this simulation 
method. 

Another object of the present invention is to provide a program 
15 which attains a two-dimensional or three-dimensional process simulation 
method which can shorten calculation time and improve the accuracy of 
calculation even in a fine region resulting from, for example, the progress 
of thin film-formation and a recording medium in which this program is 
recorded. 

20 To attain the above object, according to a first aspect of the 

present invention, there is provided a simulation method comprising a first 
step of generating plural boundary points on a string formed on the surface 
of a material, a second step of obtaining a first length of a line segment 
between the boundary points, a third step of calculating the displacement 

25 of the boundary point according to a process model, a fourth step of 

moving the boundary point by the displacement, a fifth step of obtaining a 
second length of the line segment between the boundary points after the 
boundary point is moved and a sixth step of making reference to the first 
and second lengths to add or eliminate boundary points. In this structure, 

30 the line segment after the boundary point is moved is regulated by making 
reference to the length of the stored line segment. Therefore, it is not 
required for users to design a maximum length and a minimum length to 
regulate the line segment. Further, such a problem that useless boundary 
points are increased can be avoided. This ensures that calculation time is 

35 more shortened than in a conventional method and a two-dimensional and 
three-dimensional topography simulation can be accomplished with high 



calculation accuracy. 

According to a second aspect of the present invention, there is 
provided a simulator comprising a node displacement-calculating section 
for calculating the displacement of the boundary point in accordance with 
5 a process model, a node moving section for moving the boundary point by 
the displacement, a length-calculating section for obtaining the length of 
the line segment between the boundary points before and after the 
boundary point is moved and a regulating section for making reference to 
the lengths before and after the boundary point is moved to add or 

10 eliminate the boundary point. According to this structure, the line 
segment after the boundary point is moved is regulated by making 
reference to the length of the stored line segment. Therefore, it is not 
required for users to design a maximum length and a minimum length to 
regulate the line segment. Further, such a problem that useless boundary 

15 points are increased can be avoided. This ensures that calculation time is 
more shortened than in a conventional method and a two-dimensional and 
three-dimensional topography simulation can be accomplished with high 
calculation accuracy. 

According to a third aspect of the present invention, there is 

20 provided a program which is executed by a computer, the program 
comprising a first step of obtaining a first length of a line segment 
between the boundary points generated on the surface of a material, a 
second step of calculating the displacement of the boundary point 
according to a process model, a third step of moving the boundary point by 

25 the displacement, a fourth step of obtaining a second length of the line 
segment between the boundary points after the boundary is moved and a 
fifth step of making reference to the first and second lengths to add or 
eliminate boundary points. In this structure, the line segment after the 
boundary point is moved is regulated by making reference to the length of 

30 the stored line segment. Therefore, it is not required for users to design a 
maximum length and a minimum length to regulate the line segment. 
Further, such a problem that useless boundary points are increased can be 
avoided. This ensures that calculation time is more shortened than in a 
conventional method and a two-dimensional and three-dimensional 

35 topography simulation can be accomplished with high calculation 
accuracy. 
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Other and further objects and features of the present invention will 
become obvious upon an understanding of the illustrative embodiments 
about to be described in connection with the accompanying drawing or 
will be indicated in the appended claims, and various advantages not 
5 referred to herein will occur to one skilled in the art upon employing of 
the invention in practice. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1A and FIG. IB are views for explaining the terms used in a 
10 simulation method. 

FIG. 2 to FIG. 4 are flowcharts of a conventional simulation 
method for topography processing. 

FIG. 5A and FIG. 5B are views for explaining the state of a silicon 
oxide film formed in a trench to which a simulation method is applied. 
15 FIG. 6 is a view for explaining a conventional simulation method 

in which the conventional simulation method is applied to the topography 
of a silicon oxide film formed in a trench. 

FIG. 7 is a schematic structural view of a simulator according to 
an embodiment of the present invention. 
20 FIG. 8 is a flowchart of a node-moving method in a simulation 

method according to the present invention. 

FIG. 9A to FIG. 9C are sectional views of a silicon substrate for 
explaining a simulation method for the topography of a thermal oxide film 
when a silicon substrate is thermally oxidized using a node-moving 
25 method in a simulation method according to the present invention. 

FIG. 10 is a flowchart of a simulation method according to an 
embodiment of the present invention. 

FIG. 11 is a flowchart showing the addition of boundary points, 
the flow being the internal flow of the step S26 of FIG. 10. 
30 FIG. 12 is a flowchart showing the elimination of boundary points, 

the flow being the internal flow of the step S26 of FIG. 10. 

FIG. 13 is a structural view of the data used in a simulation 
method according to an embodiment of the present invention. 

FIG. 14 is an outward view showing the relation between a 
35 recording medium and a computer according to an embodiment of the 
present invention. 
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FIG. 15A and FIG. 15B are views for explaining the state of a 
silicon nitride film (simulation example 1) formed around a gate electrode 
of a MOSFET to which film a simulation method according to an 
embodiment of the present invention is applied. 

FIG. 16, FIG. 21 and FIG. 22 are views for explaining a simulation 
method for a silicon nitride film in a simulation example 1. 

FIG. 17 to FIG. 20 are views showing changes in a data structure 
for explaining a simulation method for a silicon nitride film in the 
simulation example 1. 

FIG. 23 is a view for explaining a simulation method according to 
an embodiment of the present invention, the method being applied to the 
topography of a silicon nitride film formed in a trench of a simulation 
example 2. 

DETAILED DESCRIPTION OF THE EMBODIMENTS 

Various embodiments of the present invention will be described 
with reference to the accompanying drawings. It is to be noted that the 
same or similar reference numerals are applied to the same or similar parts 
and elements throughout the drawings, and the description of the same or 
similar parts and elements will be omitted or simplified. 

Explanations will be furnished on the assumption that the analytic 
zone is two-dimensional and the displacement components are u x (a 
displacement in the x direction) and u y (a displacement in the y direction) 
for the sake of convenient explanations. 

FIG. 7 is a schematic structural view of a simulator according to 
the present invention. A simulator 1 according to the present invention 
comprises a calculating section 2, a program recording section 3 in which 
a program to be executed in the calculating section 2 is recorded, a data 
recording section 4 which records, for example, data to be used for the 
calculation of all nodes and line segments and an input/output control 
section 5 which exchanges data with an output device 6 and an input 
device 7. The calculating section 2 comprises a node 
displacement-calculating section 11 which calculates the displacement of 
every node, a node-moving section 12 which moves a node based on the 
displacement of every node, a node-moving section 13 which sets a node 
in the initial oxidation stage, a calculating section 14 which calculates the 
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length of a line segment between nodes and a regulating section 15 which 
adds or eliminates a boundary point on the basis of the calculated length 
of a line segment to regulate the length. Further, the line segment 
length-regulating section 15 comprises a boundary point-addition section 
5 16 and a boundary point-elimination section 17. The details will be 
explained with reference to a flowchart of the following embodiment. 

(With regard to a method of moving a node) 

A topography processing method in the case of applying the 

10 node-moving method to a simulation of a thermal oxidation step will be 
explained. In the node-moving method, only the node 
displacement-calculating section 11 and the node-moving section 12 are 
used in the calculating section 2 shown in FIG. 7. 

FIG. 8 shows a flowchart of the node-moving method according to 

15 the present invention. 

First, in step SI, the displacements of all nodes No to N„-i in a 
region of a material are calculated every each of the components u x and u y 
in the node displacement-calculating section 11. The following steps S2 
to S8 are carried out in the node-moving section 12. In the step S2, the 

20 displacement components calculated in the step SI are added every 
component to the displacement components stored until the last time. 
The initial value of each displacement component is set to 0(zero). In the 
step S3, the identification number i of a node is set to 0(zero). In the 
step S4, whether the displacement {V(u 2 x0 + u 2 y0 )} of a node No is more 

25 than a displacement tolerance (l t0 ]) or not is judged. If the displacement 
of the node No is more than the displacement tolerance, the node No is 
moved in the step S5, specifically, the displacement components are added 
every each of the position coordinates of the node N 0 . Then, in the step 
S6, the displacement components (u x o, u y o) are cleared and set to O(zero). 

30 If the displacement of the node No is less than the displacement tolerance, 
the displacement component is stored without moving the node No. In 
the step S7, 1 is added to the identification number i to set the 
identification number newly to 1. In the step S8, whether or not the 
identification number 1 is smaller than a number n obtained by adding 1 to 

35 the maximum value n-1 of the identification number is judged. If smaller, 
the step S4 is carried out again. Then, this loop is repeated until the 



identification number i is n-1. When the identification number i is n in 
the step S7, the judgment made in the step S8 is "NO" and the flowchart of 
the node-moving method is finished. This method ensures that the 
accumulation of errors caused by the displacements less than the 
5 displacement tolerance can be eliminated thereby reducing errors in actual 
topography movement. 

(Simulation example by a node-moving method) 

A topography processing method in the case of applying the 

10 node-moving method to a simulation of a thermal oxidation step. In the 
simulation example, as shown in FIG. 9A, a thermal oxide film 19 is 
formed in advance on a silicon substrate 18 having a flat structure. The 
substrate 18 is thermally oxidized further. Nodes 20, 21, 22 and 23 are 
arranged only at boundary points constituting the oxide film region 19 for 

15 the sake of convenient explanations. 

First, in the step SI of FIG. 8, the displacements of all nodes 20, 
21, 22 and 23 in the oxide film region 19 are calculated every each of the 
components u x and u y . In the step S2, the displacement components 
calculated in the step SI are added every component to the displacement 

20 components initialized to 0. In the step S3, the identification number i of 
a node is set to 0. In the step S4, whether the displacement {v / (u 2 x o + 
u 2 y o)} of a node N 0 (assumed as the node 20) is more than a displacement 
tolerance (l to i) or not is judged. The displacement tolerance is, for 
example, 0.05 nm. If the displacement of the node 20 is slightly smaller 

25 than 0.05 nm due to a numerical calculation error, it is less than the 

displacement tolerance. For this, as shown in FIG. 9B, the node 20 is not 
moved and the displacement components are stored. In the step S7, 1 is 
added to the identification number i to set the identification number newly 
to 1. In the step S8, whether or not the identification number 1 is smaller 

30 than a numeral 4 obtained by adding 1 to the maximum value 3 of the 

identification number is judged. Since the judgment is true, the step S4 
is carried out again. 

In the step S4, whether the displacement of a node N\ (assumed as 
the node 21) is more than a displacement tolerance (l to i) or not is judged. 

35 Since the displacement of the node 21 is more than the displacement 
tolerance, the node N21 is moved in a step S5 as shown in FIG. 9B. 
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Actually, the displacement component is added every component of the 
position coordinates of the node N21. Then, in a step S6, the 
displacement component is cleared and set to 0. In a step S7, 1 is added 
to the identification number i to set the identification number newly to 2. 
5 In a step 8, whether the identification number 2 is smaller than 4 or not is 
judged. Since 2 is smaller than 4, the step S4 is carried out again. The 
node 22 which is a node N2 and the node 23 which is a node N3 are 
processed by the same loop as in the case of the node 21. The 
displacement of the node 23 is slightly larger than 0.05 nm due to a 

10 numerical calculation error and the node 23 is moved downward slightly 
longer than 0.05 nm from a standard plane 24. Finally, in the step S7, 
when the identification number i is 4, the judgment made in the step S8 is 
"NO" in the step S8 and the flowchart of the node-moving method is 
finished. As a result, the interface between the silicon substrate 18 and 

15 the oxide film 19 is not made to be flat. 

Next, a process according to the flowchart of FIG. 8 is performed 
again using the resulting topography of FIG. 9B as a starting topography. 
First, the displacements of all nodes 20, 21, 22 and 23 are calculated in the 
step SI. In the step S2, the displacements components calculated in the 

20 step SI are added to the displacement components which are stored or 

cleared in the last flow. In the step S3, the identification number i of the 
node is set to 0. In the step S4, whether the displacement of the node No 
(node 20) is greater than the displacement tolerance (l, 0 i) or not is judged. 
The displacement of the node 20 is a value of about 0.1 nm which is the 

25 sum of a value slightly smaller than 0.05 nm and a value of the same level. 
Since the displacement of the node 20 is more than the displacement 
tolerance, the node 20 is moved in the step S5 as shown in FIG. 9C. 
Then, in the step S6, the displacement components are cleared. In the 
step S7, 1 is added to the identification number i to set the identification 

30 number newly to 2. In the step S8, whether the identification number 2 is 
smaller than 4 or not is judged. Since the judgment is true, the step S4 is 
carried out again. 

In the step S4, whether the displacement of a node Ni (assumed as 
the node 21) is more than the displacement tolerance (l to i) or not is judged. 

35 Since the displacement of the node 21 is more than the displacement 

tolerance, the node 21 is moved in a step S5 as shown in FIG. 9C. In the 
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step S6, the displacement component is cleared. In the step S7, 1 is 
added to the identification number i to set the identification number newly 
to 2. In the step 8, whether the identification number 2 is smaller than 4 
or not is judged. Since the judgment is true, the step S4 is carried out 
5 again. The node 22 which is the node N2 and the node 23 which is the 
node N3 are processed by the same loop as in the case of the node 21. 
The displacement of the node 23 is slightly larger than 0.05 nm. The 
node 23 is moved downward slightly longer than 0.05 nm and disposed at 
a position slightly longer than 0.1 nm downward from the standard plane 

10 24. Finally, in the step S7, when the identification number i is 4, the 
judgment made in the step S8 is "NO" and the flowchart of the 
node-moving method is finished. As a result, the interface between the 
silicon substrate 18 and the oxide film 19 is made flat. 

Specifically, the displacement components of the node 20 are 

15 stored. Therefore, when the substrate is oxidized in the state shown in 
FIG. 9B, the displacement components of the node 20 are those obtained 
by adding the sum of displacement components until the last time, which 
makes it possible to be freed of the influence of numerical calculation 
errors as shown in FIG. 9C. Also, according to the present invention, 

20 such a problem that the oxide film region 19 cannot be formed even if the 
time interval is decreased to reduce the numerical calculation errors does 
not arise. 

(Embodiment) 

25 In the embodiment, a topography processing method in the case of 

applying the method to a simulation in a deposition step will be explained. 
In the embodiment, besides the node displacement calculating section 11, 
the node-moving section 12 and the node-removing section 13 in the 
initial oxidation stage, a line segment length-calculating section 14 and a 

30 boundary point addition section 16 and a boundary point elimination 

section 17 of the line segment length-regulating section 15 are newly used 
in the calculating section 2 shown in FIG. 7. 

In a topography simulation for calculating the topography of a 
semiconductor device in etching and deposition steps or the like, first the 

35 topography of the surface of a two-dimensional section is expressed 

approximately by a series of plural line segments s. Next, using a node N 
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among these line segments s as a boundary point P, the boundary point P is 
moved every small time interval in accordance with each process model to 
calculate the hourly variation of the topography. At this time, as 
topography processing technologies which become of importance, not only 
5 the technologies concerning the aforementioned movement of the node N 
but also technologies for regulating the length of the line segment s, 
specifically, technologies for adding and eliminating the node N are 
required. 

FIG. 10 shows a flowchart of a topography simulation method 

10 involving a method of regulating the length of the line segment according 
to the embodiment of the present invention. First, in a step S21, the 
number of loops for moving the boundary point P is set as the number ni oop 
of steps and 0 is set as the present number j of steps. Next, in a step S22, 
the length 1 of each of one or more line segments s constituting a string is 

15 calculated and stored in the line segment length-calculating section 14 of 
FIG. 7. In a step S23, 1 is added to the present number j of steps to set 1 
as the number of steps. In a step S24, the displacement component of the 
boundary point P is calculated in the same manner as in the step SI of FIG. 
8. In a step S25, the boundary point P is moved in the same manner as in 

20 the steps S2 to S8 of FIG. 8. In a step S26, the length of the line segment 
s is regulated. The above steps are repeated until the number j of steps 
for the calculation of the displacement is equal to ni oop . 

FIG. 11 is a flowchart showing an internal flow for adding the 
boundary point P to regulate the length r of the line segment s in the step S 

25 26 of FIG. 10. First, in a step S31, the length of the line segment s after 
the boundary point P is moved is calculated in the line segment 
length-calculating section 14 shown in FIG. 7. The following steps S32 
to S37 are carried out in the boundary point addition section 16 shown in 
FIG. 7. In the step S32, the identification number i of the line segment s 

30 is set to 1. In the step S33, whether or not the length ri of the line 

segment Si after the boundary point P is moved is more than a value axli 
obtained by multiplying the stored length h before the boundary point is 
moved by a factor a is judged. If ri > axl 1? a new boundary point P is 
added to the line segment si in the step S34. In the step S35, a length 1 is 

35 stored for each of two line segments s produced by dividing the line 

segment sj. In the step S36, 1 is added to the identification i to set the 
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number i to 2. On the other hand, if ri is smaller than axl ]; the step S36 
is carried out without carrying out the steps S34 and S35. The process 
involving the step S33 to the step S36 is repeated until the identification 
number i of the line segment s is equal to m-1. 
5 Here, the factor a is designed to be larger than 1. This makes it 

possible to attain such an object as to provide a new boundary point with 
the intention of shortening the interval between each boundary point P 
when the interval is spread so that fine topography can be obtained with 
difficulty. Also, since the line segment s is divided into two line 

10 segments by a new boundary point, the factor a is preferably about 2 to 
retain the same length as that of the original line segment. On the other 
hand, if the factor a exceeds 4, the original purpose of obtaining a detail 
topography cannot be attained. In addition, it is considered that the 
accuracy of the simulation is more improved as the factor is closer to 1 

15 from 2. 

FIG. 12 is a flowchart showing an internal flow for eliminating the 
boundary point P to regulate the length r of the line segment s in the step S 
26 of FIG. 10 similarly to the case of FIG. 11. First, in a step S41, the 
length r of the line segment s after the boundary point P is moved is 

20 calculated in the line segment length-calculating section 14 shown in FIG. 
7. The following steps S42 to S49 are carried out in the boundary point 
eliminating section 17 shown in FIG. 7. In the step S42, the 
identification number i of the line segment s is set to 1. In the step S43, 
whether or not the length ri of the line segment Si after the boundary point 

25 P is moved is less than a value |3xli obtained by multiplying the stored 

length li before the boundary point is moved by a factor |3 is judged. If ri 
s (3x1^ one of the boundary points P constituting the line segment Si is 
eliminated and a new line segment s'i is generated in the step S44. In the 
step S45, a length li of the line segment s'i is stored. Next, in the steps 

30 S46 and S47, if the eliminated boundary point P is a boundary point P 

common to the line segments Si and S2, S2 is replaced with s'i. In the step 
S48, 1 is added to the identification i to set the number i to 2. On the 
other hand, if n is larger than (3xl l5 the step S48 is carried out without 
carrying out the steps S44 and S47. Also, if the eliminated boundary 

35 point P is not a boundary point P common to the line segments si and S2, 
the step 48 is carried out without carrying out the steps S44 to S47. The 
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process involving the above step S43 to the step S48 is repeated until the 
identification number i of the line segment s is equal to m-1. 

Here, the factor |3 is designed to be smaller than 1. This makes it 
possible to narrow the interval between each boundary point P and to 
5 attain such an object as to omit calculations more than required to obtain a 
desired topography. Also, since two line segments s are made into one 
line segment s, the factor |3 is preferably about 0.5 to retain the length 
longer than that of the original line segment. On the other hand, if the 
factor f5 is less than 0.25, the original purpose of shortening the 

10 calculation time cannot be attained. In addition, it is considered that the 
calculation time is more shortened as the factor is closer to 1 from 0.5. 
As the factor is closer to 0.5 on the contrary, the accuracy of the 
simulation is more improved. It is to be noted that if |3 is designed to be 
an inverse of a, the variables can be decreased with the result that the 

15 factors are reduced to one in number to attain the optimization of the 
simulation with ease. 

FIG. 13 is a view of a data structure used in a simulation method 
according to the embodiment of the present invention. The data structure 
comprises a boundary point identification number data zone 41 to which 

20 the identification numbers of all boundary points and line segments can be 
input, a position coordinate data zone 42 to which the coordinates of a 
boundary point can be input corresponding to an identification number, a 
boundary point displacement data zone 43 to which the displacement of a 
boundary point can be input every coordinate component corresponding to 

25 an identification number, a before-move line segment length data zone 44 
to which the length of a line segment before a boundary point is moved 
can be input corresponding to an identification number, an after-move line 
segment length data zone 45 to which the length of a line segment, after a 
boundary point is moved and before the length is regulated, according to 

30 an identification number and a line segment length data zone 46 to which 
the length of a line segment after the length is regulated can be input 
corresponding to an identification number. 

The series of processings shown by a flowchart in FIG. 8 and FIG. 
10 to FIG. 12 is recorded in a computer-readable recording medium as a 

35 program to be executed by a computer. The data structure having data 
operated by a computer as shown in FIG. 13 is also recorded in the 
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computer-readable recording medium. This program is stored in the 
program-recording section 3 of FIG. 7 and the data structure is stored in 
the data recording section 4. Here, the recording medium includes 
mediums capable of recording a program such as a semiconductor memory, 
5 magnetic disk, optical disk and magnetic tape. 

FIG. 14 is an outward view showing one example of a computer 
system which reads a program and the like stored in a recording medium 
according to the embodiment of the present invention to carry out a 
topography simulation according to the procedures written therein. A 

10 floppy disk drive 31 and a CD-ROM drive 32 are disposed on the front of 
the body of a computer system 30. A floppy disk 33 or a CD-ROM 34 
are inserted from the entrance of each drive and a prescribed read 
operation is carried out to thereby install programs stored in these 
recording media in the system. Further, the program can be transferred 

15 from other computers through a transmission medium such as internet. 

Also, for example, a ROM 35 as a semiconductor memory used for a game 
pack and a cassette tape 36 as a magnetic tape may be used by connecting 
specified drives. 

20 (Simulation example 1 of the embodiment) 

Explanations will be furnished as to a topography processing 
method in the case of applying the embodiment to a simulation of a step of 
deposition of a silicon nitride film for forming a side wall of a field effect 
transistor (MOSFET) having a metal oxide semiconductor structure. In 

25 the simulation example, as shown in FIG. 15A, a thermal oxide film 19 to 
be a gate oxide film is disposed on a flat silicon substrate 18 and a 
polysilicon 25 to be a gate electrode is disposed on the film 19. A square 
26 is a simulation region which is a subject of calculation this time. As 
shown in FIG. 15B, a silicon nitride film 27 is formed on the exposed 

30 surface of each of the substrate 18, film 19 and polysilicon 25 shown in 
FIG. 15A. This filming step is simulated hereinbelow to obtain the 
topography of the film 27 to be generated in the region 26. 

Flowcharts shown in FIG. 8, FIG. 10 and FIG. 11 are used for the 
simulation. The chart of FIG. 10 is used as a main routine and the charts 

35 of FIG. 8 and FIG. 11 are used as subroutines. First, in a step S21 of FIG. 
10, the number ni oop of steps for the calculation of a displacement is set to 
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3 and the present number j of steps is set to 0. 

Next, in FIG. 16, a string which connects the solid line through □ 
is generated. Here, the string means a plane to be observed for the 
purpose of a simulation. The position of the solid line is made to accord 
5 to the surface of the polysilicon 25. □ represents a boundary point P 
and identification numbers 0 to 6 are given to the boundary points P 
respectively. The boundary point P is arranged at higher densities on 
solid lines considered to be significantly changed in topography. In the 
simulation example 1, many boundary points P are arranged in the vicinity 

10 of the corners of the polysilicon. A solid line between each boundary 
point P is a line segment s. Identification numbers 0 to 6 are also given 
to the line segments s respectively. In the step S22, the respective 
lengths h to U of six line segments s are calculated and stored. At this 
time, the data is stored in the data structure as shown in FIG. 17. In the 

15 zone 41, identification numbers 0 to 6 are placed. In the zone 42, 
position coordinates (x, y) of every boundary point P are input every 
component corresponding to each number. In the zone 44, the lengths li 
to 1 6 are input corresponding to each number of the line segments s. In 
the zone 43, since the displacements are accumulated, all displacement 

20 components in the zone was reset to 0. 

In the step S23, 1 is added to the present number j of steps to set 
the number of steps to 1. In the step S24, displacement components Ux 
and Uy of each of the boundary points Pi to P6 are calculated. As shown 
in FIG. 18, the calculated displacement component is added to the already 

25 reset value 0 every number and every component in the zone 43. In the 

step S25, the boundary point P is moved in the same manner as in the steps 
S2 to S8 shown in FIG. 8. In this simulation example 1, because each 
displacement of the boundary points Pi to P6 is larger than the 
displacement tolerance (l to i), these boundary points Pi to P6 are moved as 

30 shown by the arrow in FIG. 16 to obtain each A mark. In the data 
structure, as shown in FIG. 18, the displacement component is added to 
the position coordinate component every number and every component. 
To state in more detail, in the case of the position coordinate x of the 
number 0, xo' is found by calculation from the equation: Xo + u x o = Xo' and 

35 the obtained x 0 ' is stored in the zone of the position coordinate x of the 

number 0 as shown in FIG. 19. The above process is based on the step S5 
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of FIG. 8. In succession, in the step S6, the zone 43 is cleared and set to 
0 as shown in FIG. 19. 

In the step S26 shown in FIG. 10, the length of the line segment s 
is regulated. For the regulation of the length, the flowchart of FIG. 11 
5 used to add a boundary point is used. In the step S31, the respective 
lengths ri to re of the line segments si to S6 after the boundary point P is 
moved are calculated. Then, the lengths ri to t6 are arranged in the zone 
45 of FIG. 19. In the step S32, the identification number i of the line 
segment s is set to 1. In the step S33, whether or not the length ri is 

10 more than a value ctxl] obtained by multiplying the length li before the 
boundary point is moved by a factor a is judged. In this simulation 
example 1, a is assumed as 1.8. Since the equation ri > axlj is not valid, 
next the step S36 is carried out, so that the identification number i is set to 
2. In the step S37, the flow is returned to the step S33. The process 

15 involving the step S33 to the step S36 is repeated until the identification 
number i of the line segment s is equal to 6. 

In this simulation example 1, when the identification number i is 3 
and 4, ri > oixli. For example, when the number i is 3, the judgment in 
step S33 is YES. So, a new boundary point P 3 is added to the line 

20 segment S3 in the step S34 and the lengths I3' and I4' of each of two line 
segments s generated by dividing the line segment Si are stored in the step 
S35. In the data structure, as shown in FIG. 20, the coordinates x 7 ' and 
y 7 ' of the new boundary point P 3 are generated and the generated 
coordinates are inserted in the zone of the coordinates of the number 3. 

25 When these coordinates are inserted, the zones of the position coordinates 
of the series of numbers successive to and including the number 3 are 
displacemented to the zones of the series of numbers successive to and 
including the number 4. The lengths I3' and U' are arranged in the zone 
46. The line segment S4 is subjected to the same process to obtain the 

30 dotted line of FIG. 16 and the data structure of FIG. 20. Finally, the data 
of the zone 46 is overwritten on the zone 44 to make the state of the data 
structure of FIG. 20 equal to that of the data structure of FIG. 17. 

The step S26 of FIG. 10 is finished and in the step S27, whether or 
not the present number of steps, namely, 1 reaches 3 which is the number 

35 nioop of steps is judged. Since the present number does not reach the 
number n] 00p , the flow is returned to the step S23 to set the number j of 



steps to 2. Thereafter, the steps S24, S25 and S26 are carried out in the 
same manner as in the case where j = 1 to obtain the solid line segment 
connected through □ and O as shown in FIG. 21. Finally, j is set to 3 
and the same process is carried out in the same manner as in the case 
5 where j = 1 to obtain the dotted line segment connected through A as 
shown in FIG. 22. This dotted line segment shows the topography of the 
film 27 shown in FIG. 15B. It is to be noted that the change of j signifies 
a transition of time and the string obtained by connecting a line segment 
obtained every j shows the topography obtained at the time corresponding 
10 to j. By applying the present invention, the length of a line segment after 
a boundary point is moved can be most properly regulated and the 
topography of the nitride film 27 at the corner portions of the gate 
polysilicon can be simulated with high accuracy. 

15 (Simulation example 2 of the embodiment) 

Explanations will be given as to a topography processing method 
in the case of applying the embodiment to a simulation of a step of 
depositing a silicon oxide film on a groove formed on a silicon substrate. 
In the simulation example 2, as shown in FIG. 5A, a groove 59 is formed 

20 on a silicon substrate 52. A square 60 represents a simulation region 
subject to calculation this time. As shown in FIG. 5B, a silicon oxide 
film 61 is formed on the surface of the silicon substrate 52 including the 
internal surface of the groove 59. Hereinafter, this filming step is 
simulated to determine the topography of the film 61 generated in the 

25 region 60. 

Flowcharts used for the simulation are FIG. 8, FIG. 10, FIG. 11 
and FIG. 12. The chart of FIG. 10 is used for a main routine and the 
flowcharts of FIG. 8, FIG. 11 and FIG. 12 are used for a subroutine. 
Because, in this example 2, the flowcharts of FIG. 11 and FIG. 12 are used 

30 to regulate a length, these flowcharts of FIG. 11 and FIG. 12 may be 

united to form a new flowchart, which may also be used in this example 2. 
Specifically, the steps S31, S32, S36 and S37 of FIG. 11 and the steps S41, 
S42, S48 and S49 are unified respectively and the obtained steps are 
defined as steps S31, S32, S36 and S37. Then, the cases of judgment in 

35 the step S33 of FIG. 1 and in the step S43 of FIG. 12 are united into three 
cases of judgment. A first case is the case where the judgment is YES in 
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the step S33. In this case, the steps S34, S35 and S36 are carried out in 
this order. A second case is the case where the judgment is YES in the 
step S43. In this case, the steps S44, S45, S46, S47 and S36 are carried 
out in this order. Finally, a third case is the case where the judgment is 
5 NO in any of the steps S33 and S43. In this case, the step S36 is carried 
out eventually. 

First, in the step S21 of FIG. 10, the number ni oop of steps for the 
calculation of a displacement is set to 1 and the present number j of steps 
is set to 0. 

10 Next, a string connecting the solid line of FIG. 23 through □ is 

generated. The position of the solid line is made to accord to the surface 
of the substrate 52. □ represents a boundary point P. Identification 
numbers 0 to 16 are respectively given in order to each boundary point P. 
Each boundary point P is arranged on the solid line in closer contact at the 

15 side surface of the groove 59 and its vicinity where the topography is 

considered to be significantly changed. On the contrary, each topography 
between the boundary points P0 and PI and between the boundary points 
P15 and P16 is considered to be able to be estimated from each 
displacement between the boundary points P0 and PI and between the 

20 boundary points P15 and P16 respectively and no node is arranged 

between the boundary points P0 and PI and between the boundary points 
P15 and PI 6, whereby it is expected to shorten calculation time. The 
solid line between each boundary point P is a line segment s. 
Identification numbers 1 to 16 are also given to each line segment s. In 

25 the step S22, the respective lengths li to 1 5 6 of 16 line segments s are 
calculated and stored. 

In the step S23, 1 is added to the present number j of steps to set 
the number of steps to 1. In the step S24, the displacement components 
Ux and Uy of each of the boundary points Pi to Pi6 are calculated. In the 

30 step S25, the boundary point P is moved in the same manner as in the steps 
S2 to S8 shown in FIG. 8. In this example, each displacement of the 
boundary points P] to Pie is larger than the displacement tolerance (l to i), 
these boundary points Pi to Pi 6 are moved as shown by the arrow in FIG. 
23 to obtain each A mark. 

35 In the step S26 shown in FIG. 10, the length of the line segment s 

is regulated. When the length is regulated, the above flowchart obtained 
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by uniting the flowchart of FIG. 11 with the flowchart of FIG. 12 is used. 
In the step S31, the respective lengths ri to rie of the line segments Si to 
S]6 after the boundary point P is moved are calculated. In the step S32, 
the identification number i of the line segment s is set to 1. In the steps 
5 S33 and S43, whether the length ri is more than axlj (the first case) or 
less than |3xli (the second case) or neither more than axlj nor less than 
|3xl! (the third case) is judged. In this example, a is set to 1.2 and |3 is 
set to 0.5 and hence, ri fits the third case. Therefore the step S36 is then 
carried out and the identification number i is set to 2. In the step S37, 

10 the flow is returned to the step S33. The process involving the above 
steps S33 to S36 is repeated until the identification number i of the line 
segment is equal to 16. 

In this simulation example 2, the first case where r is larger than 
axl is the case where the identification number i is 13 and 14. For 

15 example, in the case where i is 13, the judgment is YES in the step S33, a 
new boundary point 28 is added to the line segment S13 in the step S34 as 
shown in FIG. 23 and the lengths of each of two line segments s obtained 
by dividing the line segment S13 are stored. 

In this simulation example 2, the second case where r is smaller 

20 than (3x1 is the case where the identification number i is 3 and 4. For 
example, when the number i is 3, the judgment is YES in the step S43. 
Then, in the step S44, the boundary point P3 constituting the line segment 
S3 is eliminated and a new line segment s'3 is generated. In the step S45, 
the length I3 of the line segment s'3 is stored. Next, in the steps S46 and 

25 S47, if the eliminated boundary point P3 is a boundary point common to 
the line segments S3 and s 4 , s 4 is replaced with s' 3 . 

The step S26 of FIG. 10 is finished and in the step S27, whether or 
not the present number of steps, namely, 1 reaches 1 which is the number 
nioop of steps is judged. Since the present number reaches the number 

30 nioop, the flow of FIG. 10 is finished. Finally, a line segment shown by 

the dotted line connected through A and O is obtained as shown in FIG. 
23. This line segment shown by the dotted line shows the topography of 
the film 61 of FIG. 5B which is intended to be found. 

As shown in FIG. 23, the boundary points 62 and 64 shown in FIG. 

35 6 are not generated in FIG. 23. This makes it possible to shorten the 
calculation time. On the other hand, a boundary point which is not 
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generated on the line segment S13 in FIG. 6 can be generated as the 
boundary point 28 in FIG. 23. This enables the accuracy of the 
topography to be improved. Also, a reduction in processing hours and an 
improvement in accuracy can be attained at the same time. 
5 In the embodiment, the length of a stored line segment is reset for 

a new line segment generated by the addition or elimination of a boundary 
point. However, it is possible to improve the accuracy of calculation by 
calculating the length of the generated new line segment and setting the 
calculated value to regulate the length of the line segment. Also, in the 
10 case of a three-dimensional simulation, the area of a polygon may be used 
in place of the length of a line segment. Moreover, the present invention 
can be applied to topography simulations in steps other than a deposition 
step. 

Various modifications will become possible for those skilled in the 
15 art after receiving the teachings of the present disclosure without 
departing from the scope thereof. 
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